Evidence of several dipolar quasi-invariants in Liquid Crystals 
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In a closed quantum system of N coupled spins with magnetic quantum number I, there are 
' about (21 + 1) N constants of motion. However, the possibility of observing such quasi-invariant 

(QI) states in solid-like spin systems in Nuclear Magnetic Resonance (NMR) is not a strictly exact 
prediction. The aim of this work is to provide experimental evidence of several QI, in the proton 
1 NMR of small spin clusters, besides those already known Zeeman, and dipolar orders (strong and 

weak). We explore the spin states prepared with the Jeener-Broekaert pulse sequence by analyzing 
the time-domain signals yielded by this sequence as a function of the preparation times, in a variety 
of dipolar networks. We observe that the signals can be explained with two dipolar QIs only within 
a range of short preparation times. At longer times the time-domain signals have an echo-like 
behaviour. We study their multiple quantum coherence content on a basis orthogonal to the z-basis 
and see that such states involve a significant number of correlated spins. Then we show that the 
whole preparation time-scale can only be reconstructed by assuming the occurrence of multiple QI 
CZ) , which we isolate experimentally. 

"o . 

I. INTRODUCTION 

Quasi-equilibrium states (QE) found in Nuclear Magnetic Resonance (NMR) of strongly interacting nuclear spin 
systems are quantum states which do not evolve under the system Hamiltonian, thus they can be represented by a 
reduced spin density operator diagonal in blocks in the eigenbasis of the spin-environment Hamiltonian P, Q ■ These 
states only evolve due to spin-lattice relaxation, over a time scale much longer than the one of the build-up of the quasi- 
equilibrium. An experimental NMR procedure^to prepare and detect dipolar QE states in high-field solid-state NMR 
is the Jeener Broekaert pulse sequence (JB) [3j. Briefly, this technique consists of the phase-shifted radiofrequency 
Oh| (rf) pulses: 90^ — r — 45j, — t e — Ab y — t. The first pulse creates single quantum coherences in the spin system, which 
evolve in the rotating frame mainly under the dipole spin-spin Hamiltonian during r. Along this period, multi-spin 
,_H ' single-quantum coherences can develop, and the second 45 y pulse transforms part of the coherences just created into 
, multi-spin order 0,0, Two important processes occur during the evolution period t e , along different timescales: 
decoherence and relaxation. Finally, the third pulse converts QE states into observable single quantum coherence. 

JB experiments are frequently used to transfer the Zeeman order into dipolar order, which, in a variety of samples 
is the only observable dipolar QE state. It is an experimental fact that for short preparation times r <C 1/oJd m 
the JB sequence, the dipolar signal observed after the third pulse is proportional to the time derivative of the FID 
signal [l2j], which indicates that the prepared state is very similar to the secular (high field) dipolar energy. Two 
dipolar QE states, were instead observed both on hydrated salts @, 01, an d on liquid crystals (LC) 0,H|. In the 
former, the occurrence of 'intra-pair' and 'inter-pair' QE states was associated with the distribution of the proton 
. . ' pairs of the water molecules in the lattice, conforming a system of weakly interacting spin pairs. The spin system 
of LCs is composed by the few proton spins within each molecule, which display a hierarchy of dipolar couplings. 
This kind of dipolar network is different from hydrated salts and a model of spin pairs does not strictly proceed, 
however, the dipolar Hamiltonian can in fact be partitioned into two mutually commuting and orthogonal parts, the 
'strong' %s and the 'weak' Hw terms [10]. Since these operators are chosen so that they commute with the total 
spin Hamiltonian, they are nearly constants of the motion (for evolution times t e much greater than 1/wn [HI]). For 
this reason, they are called quasi-invariants (QI). 

In early works on dipolar order relaxation in LCs, the experiment was interpreted in terms of a single quasi- 
invariant, ps, of an ensemble of isolated representative stron gly coupled spin pairs however, at that time, 
the dependence of relaxation times on r was not understood |14j . Later work on LC showed that the state pw 
prepared by setting the preparation time to the one at which the ps order vanishes, also presents properties of quasi- 
equilibrium. In fact, it relaxes exponentially towards equilibrium, with a characteristic time which is also of the order 
T\d 0]- Besides, experiments where the multiple quantum coherences of the states prepared with the JB sequence 
are encoded in orthogonal, X basis, carried out in nematic 5CB [ll|, confirmed that the tensor structure of the iS 
state corresponds to two-spin dipolar order while the W state shows high order quantum coherences, all of which 
relax towards equilibrium with the same decay constant Tiyy. This indicates that for short preparation times the 
two-spin order prevails while higher order spin correlations are growing. This second QI, produces discernible signals 



2 



coresponding to multi-spin order, created within a time window r where the dipolar order vanishes, similarly to the 
hydrated crystal studied in reference Q ■ 

Recently, it was shown that quantum decoherence can provide an efficient mechanism by which the spin system 
attains a diagonal state in the basis of the spin-environment Hamiltonian, over a time scale which is intermediate 
between those governed by its own interactions, and thermalization ruled by thermal fluctuations of the environment 
[3, SB]. All these results support the assumption that in LC's, for times t e greater than the decoherence time scale, 
the density operator can be written in the form [Toj 

p D = l-PzH z -0s(r)'Hs-p W (T)n w . (1) 

In fact, the experimental behaviour of the dipolar signal of 5CB can be described with Eq.([T|) for preparation times 
t < 80/xs Q . Relaxation brings the system to thermal equilibrium with the whole surrounding world over a much longer 
time scale, t e 3> t, through processes involving energy exchange between the spins and the lattice. The characteristic 
lifetime T\r> of QE states (dipolar order relaxation time) is comparable to the common Zeeman relaxation time. 

Theoretically, in a cluster of N dipole interacting spins 1/2, as the protons of a typical LC molecule, there are at 
least 2 N constants of motion, the exact number depends on the degeneration of the dipolar Hamiltonian 2] . That 
is, such a number of spin operators are needed to span the commutative space (or diagonal in blocks space) of the 
dipolar Hamiltonian. The results obtained so far: the occurrence of two dipolar QI, are consistent with this idea, but 
accordingly, it should be possible to prepare new quasi-equilibrium states from the initial Zeeman order, which should 
be observable for longer preparation times in the JB experiment, as long as the JB signal is detectable. 

The aim of this work is to explore the spin states prepared for long preparation times in the JB experiment, looking 
for experimental evidence of multiple QIs that can be expected for dipole coupled spin clusters. The analysis is based 
on an exhaustive study of the time-domain signals in a variety of dipolar networks. 

Providing insight on the physics nature of the quasi-invariants of a spin cluster can be useful both for applications 
as for basic research. QE states are relevant observables of the spin system, providing relaxation parameters useful 
to study molecular motion in LC mesophases through their dependence on temperature and magnetic field 0, 17. 18 1 
They have been proposed as initial states for the excitation of multiple quantum coherences (MQ) in MQ-NMR 19 21 



and used in spin counting experiments I5 t Il5| . Besides, these states have been proposed as alternative to the Zeeman 
order in Magnetic Resonance Imaging [22j . Implementation of noiseless quantum memories and multi-spin quantum 
register relies on the possibility of manipulating multi-spin correlated entities which are unperturbed by decoherence 
processes [2314251 [28j | . From a basic viewpoint , the physics of systems with few degrees of freedom coupled to a quantum 
environment attracts today's attention of a widespread community because of the potentiality for applications such 
as quantum devices and quantum information processing (26l - |28j and, significantly, also because these systems arc 
testbeds for studying fundamental aspects as irreversibility, equilibration and thermalization 0, [29l |30| . 

Section ILTl contains the definition of the QI operators and their relation with the NMR signals in the JB experiment. 
In Section Mil we present an experimental survey of the manifestation of multiple quasiinvariants and propose a 
method to isolate them in samples with different geometries. The nature of the different QI is examined by means of 
relaxation experiments and by studyng the multiple quantum coherence content on the x-basis. 

II. QUASI INVARIANTS 

The spin Hamiltonian of resonant nuclei in ordered systems like solids and liquid crystals has an important contri- 
bution from the dipole-dipole coupling energy besides the Zeeman term, 

n s = n z + n D . (2) 

where the Zeeman energy in units of h is Hz = — w I z , with lj q the Larmor frequency. In the rotating frame description 
[3l| . the time evolution of any spin state is mainly driven by the dipolar Hamiltonian and often only by its secular 
part (high field approximation), 



n = ^E^ T 20- ( 3 ) 



i<3 

In liquid crystals the secular part of the dipolar coupling between nuclei i and j is 




3 



with rij the distance between spins, Oij the angle between the internuclear vector and the magnetic field and the angle 
brackets stand for an average over internal and reorientational molecular motions. The residual dipolar coupling in LC 
only involves protons within a molecule, then, the sum of Eq.Q runs over protons within a molecule. is the zero 
component of a normalized irreducible spherical tensor of rank two, which in terms of the spin angular momentum 
operators is [H, [33[ 



20 ~ V6 



Then, an arbitrary initial state po of a closed system evolves, according to the Liouville equation, like 

pit) = e- m °^ po e***. 

However, the interaction of a finite quantum system with the environment within the essentially adiabatic regime 
[l[ (or quasi- isolated regime [2!|) attenuates the off-diagonal elements of the spin state (in the basis of the spin- 
lattice interaction Hamiltonian) while preserving the diagonal in blocks. We say that a spin system is in a state of 
quasi-equilibrium when its density operator in the rotating frame has a diagonal form as 

i-^/WqJ, (4) 

where Qk are zero-trace operators which commute with the spin Hamiltonian in the rotating frame, [Qi,HP D ] = and 
satisfy the orthogonality relations 

1r{QkQk'} = 4,fc'Tr{Q fe Q fe }, 

the coefficients /3fc(r) are 

P k (T) = -Tr{Q k Af p qe (T,t)}/7r{Ql}, 

and J\f — Tr{l} is the normalization factor. 

In the JB experiment, such a state of internal quasi-equilibrium is reached when a time after the second pulse t e 
longer than decoherence has elapsed. In this way, the coefficients /3k can be interpreted as a generalization of the 



inverse spin temperatures 33j. Since the QE states can be written in terms of the constants of motion {Qk} these 
states only evolve in the long timescale of spin-lattice relaxation, so we call them quasiinvariants (QI) and their 
expectation values satisfy /3fc(r) = ( Qk )/Tr{Q^}. 

The NMR signal detected in phase with the read pulse, compatible with Eq.Q is 

S(t) = ( I y (t) ) = Tr{U(t) P r Pqe Pt [7(t)t I y } 

= Ei (3 t (T)1r{U(t)P r Q^UHt)I y } (5) 



= Eift(r)^(t), 

where P r represents the read pulse and U(t) — exp{— iW^t} is the evolution operator during the observation period. 
Within this view, coefficients /^(t) can be interpreted as the weight with which operator Qi contributes to the overall 
NMR signal for a given preparation time, and !Fi{i) as the contribution of the i-th quasiinvariant to the dipolar signal. 
In fact, the Zeeman QI has no projection on the signals detected on the "dipolar channel" (channel y in Eq.©) after 
the JB sequence, which is equivalent to considering /3z—0. It is worth to mention that S(t) is symmetric with respect 
to times r and t. This property, of the JB sequence. In the case that the only QI is the total dipolar energy 1-1%, as 
is the case of systems with regularly distributed nuclear spins, the dipolar signal of Eq.Q is proportional to the time 
derivative of the FID signal [33| • 

The dipolar signals for r < 80ps in both gypsum single crystal and nematic 4'-peniyZ-4-biphenyl-carbonitrile (5CB) 
were accounted for in terms of two dipolar QI. The constants of motion in gypsum are the dipolar interaction energy 
of the spin pairs (T~L°d^[) and the interaction energy of each spin with all but its pair spin, ( T-L° D z ) [H, l37j . while 
the QI used in 5CB are the strong ( Us ) and weak ( "Hyv ) parts of the dipolar interaction [3, [13, EE] ■ 
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III. EXPERIMENTAL 

The interest herein is to explore the manifestation of multiple QIs after the Jeener-Broekaert rf pulse sequence. 
With this aim we analyze the JB signals as a function of the preparation time in different kinds of dipolar networks, 
with particular attention on long preparation times. Experiments at 20MHz were carried on in a Bruker Minispec 
mq20 and those at 60MHz in a homemade pulsed NMR spectrometer based on a Varian EM360 magnet [H], IH| • 

As shown in FigflJ the response of the different compounds to variations of the preparation time is diverse, however 
they all share the characteristic of being symmetric with respect to r and t. In powder adamantane (FigQJ,), the 
signal is proportional to the time derivative of the FID and keeps this shape for every preparation time r, which is 
consistent with the occurrence of only one dipolar constant of motion: the dipolar energy, ( H° D ) [H, [36| . 

On the contrary, in the other compunds shown in FigQ] [(b) Gypsum and (c) N-(4-Methoxybenzylidene)-4- 
butylaniline (MBBA)], the dipolar signal shapes change drastically with r, being proportional to the time derivative 
of the FID only for short r (more precisely, the proportionality holds while r < to, where to is the time at which the 
time derivative of the FID first crosses through zero). 

However, a new behaviour of the dipolar signals arises clearly in LC for longer preparation times. The observed 
signals show an "echo" like behaviour, that is, they have a peaked shape whose maximum occurs at a time t from the 
read pulse equal to the preparation time r. Such behaviour is not observed in gypsum. 

In order to highlight this feature, Figs. [2]and[3]show the time t max at which the dipolar signal attains its maximum 
amplitude in the detection period, as a function of the preparation time r, in different compounds. Preparation 
times were varied within the range where the S/N ratio is greater than 1% on each compound. It can be seen that 
the experimental t max (r) (solid circles in Fig. [5]) is very similar in all the cyanobiphenyl samples, 4'-n-4-biphenyl- 
carbonitrile where n stands for pentyl (5CB), hexyl (6CB) and octyl (8CB). Data from nematic MBBA at the bottom 
of Fig. [3] behave just as the other nematic LC samples. This behaviour is drastically different in the gypsum single 
crystal (top of Fig. [3]) since t max {T) is bounded in gypsum while grows linearly for long t in all the nematic LC 
samples. As seen in this figure, the whole observable range on gypsum is rather smaller than that of the LCs, and the 
remarkable difference is that LCs develop the slope-one behaviour at preparation times near one fourth of their whole 
observable ranges while gypsum never does. It is also worth to note that the behaviour of t max (r) is also different 
in two orientations of the gypum single crystal with respect to the external magnetic field. The wavy curve of solid 
circles in Fig. [3](top) corresponds to the orentation B | [010] while the open squares correspond to an orientation 
that makes the protons at the water molecules equivalent. 

It has been already shown that two independent QI (S and W) can be prepared in nematic LC and some hydrated 
salts [t| |37j, with the property that W emerges when the weight of S becomes zero. In fact, the set of dipolar 
signals in 5CB prepared with r < 80/^s were satisfactorily reproduced by using such two QI in Eq.([S|) 9]. However, 
this procedure does not describe the signal behaviour for r > 80/xs. Now, we propose a procedure to extract the 
quasiinvariants from the experimental data, under the sole assumption that the states prepared in the experiment 
are described by Eq.Q. It might be expected that new QI, if observable, will arise sequentially as the preparation 
time increases and that in a favourable condition, a subset of preparation times {r^} exists at which one QI (the i-th) 
prevails over the others. This amounts to proposing that it is experimentally possible to isolate a subset of QI. The 
validity of such working hypothesis will be tested along this work. 

The measured data sets are composed by the signals obtained for each of the m chosen values of the preparation 
time r = T m . Let us call M Tm (t) to the acquired time domain signals after the JB sequence, which we arrange as 
the rows of the data array, and call Aft(r) to the columns of such array (pseudo-signals). Any cut A/j^r) at t = t\ is 
found to be strictly symmetric with the signal M Tl (t) prepared with t\ =t\. 

With the aim of deriving the functions J°i(t) from the experimental data we followed these steps: (i) Recognize the 
preparation time r = Ti which makes A4 Tl (t) proportional to the time derivative of the FID signal and has maximum 
amplitude. This signal attains its maximum at t = t\ = n. Based on former evidence that the state prepared with a 
time r = n is a QI, we assume that the experimental function A4 Tl (t) can be identified with the contribution of the 
first QI to the signal, T\{t) (see dotted line in Fig. 2]). (««) Also select the experimental pseudo-signal A/j^r), which 
corresponds to a cut on the data array at an observation time t = t\, and identify it with the "weight" /3i(r) of the 
first QI. If the only dipolar QI is the dipolar energy, as is the case of adamantane, the dipolar signal of Fig. [TJa) can 
be adequately reconstructed by calculating the 2D function 

Si(r,i) =Px(j)T x (t) = N tl {r)M Tl {t) 

However, this form of Si (r,t) does not describe the experimental signals of the other studied samples, then, (iii) In 
order to select the contribution of the next QI, ^(t), we analyze the difference between the measured data and the 
contribution of the first QI to the signal 



C 2 (r, t) =M T {t)- Si (r,t). 
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This data set attains its maximum value at a time which we call t% — t%. Then we choose Tiif) = C 2 (r 2 ,£) and 
/3 2 (r) = C 2 (t, t 2 ), the solid curve in Fig. S) It is worth to mention that A/t 1 (r 2 ) — 0, and therefore r 2 is the time at 
which the the first dipolar QI has negligible contribution to the observed signal (This feature agrees with previous 
works [§, Qllzj]). 

This election of T\ and T<i for reconstructing the NMR signal S 2 (r, t) yields the curves in open triangles in Figs. 
[2] and [3] (blue triangles on the online version) . It allows for an excellent reconstruction of the experimental curve 
in gypsum. In fact, as seen in the top of FigJ3] both curves (solid circles and open triangles) coincide in almost all 
the timescale. The dipolar signal shape had already been accounted for within a restricted interval of preparation 
times by Dumont et al. 37J by assuming the occurrence of two QI. Now we learn that these two QI can give a good 
description within an extended timescale. 

On the contrary in LC samples, the triangle-curves only reproduce the experimental behaviour of t max (r) for times 
r, t < 120/iS in 5CB, 6CB 8CB and MBBA, that is, when the preparation times are restricted roughly to within the first 
half period of the dipolar coupling. This failure in describing t max (r) for longer times with only two quasiinvariants 
led us to continue with the procedure in order to find a new curve F 3 . Again we analyze the residue 



C 3 (T,t)=M T (t)-S 2 (T,t) 



to find the time £3 = t 3 at which this data set attains its maximum amplitude and define F 3 (t) = C 3 {r 3 ^t) (see 
dashed curve in Fig. [3J and /^(r) = C 3 (r, t 3 ). The times t\, t 2 and i 3 corresponding to the different compounds are 
summarized in the following Table. 



Compound 


Temp (K) 


ti(fJS) 


t 2 {fjs) 


t 3 (M s ) 


5CB 


302 


30 


70 


134 


5CB 


297 


26 


60 


106 


6CB 


297 


36 


85 


160 


8CB 


311 


29 


74 


140 


gypsum 


311 


10 


28 




MBBA 


311 


22 


65 


120 



With this selection it was possible to improve the reconstructed curve for t max (r). Figures [5] (a, b,c) and[3](b) show 
that involving a third quasiinvariant yields the open circles curve which gives a noticeably better reconstruction of 
the experimental curves. FigJS] shows a detail of t max {r) in MBBA, where the three-QI curve is definitely better than 
the two-QI curve especially within the interval 115 /j,s < r <136 /us; however it still does not account for the echo-like 
behaviour for longer preparation times. 

In order to test if J~ 3 (t) does in fact behave as the signal of a QI, we analyze the signal amplitude attenuation 
as a function of the evolution time t e as in a regular spin-lattice relaxation experiment. Figure [6] shows amplitude 
attenuation of J~ 3 {t) in 5CB at 302 K, for the preparation times shown in Table 1: r = 30/is, circles; r = 70/is, open 
squares and r = 134/is, full squares. The three curves can be adequately fitted with two exponential decays, the 
common short time decay (the same for all) is 4 ms while the longer characteristic times are Tbi = 67 ms (circles), 
Td2 = 56 ms (open squares) and To 3 = 39 ms (full squares). The fast decay is just a witness of the attenuation 
of higher order coherences, in fact, its characteristic time coincides with the timescale of irreversible decoherence as 
measured in ref. [l6| . On the other hand, the longer decay times are all of the order of the dipolar relaxation time 
(67 ms in this case). 

In ref. (l5| it was shown that the quasi equilibrium state which is proportional to T-Ls in 5CB (prepared by setting 
ftw( T i) — 0) is a two-spin correlated state since it involves at most two-spin tensors, while states proportional to "Hyy 
(prepared by setting fls( T i) — 0) involve correlations of more-than-two spin tensors. With the aim of studying the 
correlated nature of the states prepared at different t's in 5CB, we analyze the multiple quantum coherence (MQc) 
content on the X-basis using the same pulse sequence as in refs. [H, HE)]. In this experiment, rotating the state around 
an axis orthogonal to Z allows encoding MQc which reflect the number of multiply connected spins in the prepared 
states. The experiment begins by preparing the state with the phase shifted JB pulse pair, then, a waiting time t w is 
followed by two pulses 9O(0 +7r / 2 )- e - 90( y ) which encode the coherence numbers of the quantum state at time t w in 
the X-basis, when varying <fr systematically in succesive experiments. Fourier transformation of the signal amplitude 
with respect to <f> yields an X-basis coherence spectrum. Our experiment, conducted on 5CB at 297 K and 60 MHz, 
was set to encode up to 8-quantum coherences on the X-basis, however, coherences higher than 4 in 5CB were seen to 
fall below the noise level. The MQc content varies with the preparation time r, as shown in Fig. [3(a). As expected, 
zero and double quantum coherences are dominant for < r < 50/xs (notice that ti in 5CB change with temperature, 
as evident in the first two rows of Table 1). The fact that the quotient of their amplitudes is ca. 1.5 indicates that 
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the states prepared within this interval can be described with only one QI (or dipolar constant of motion) Q , which 
is proportional to T 2 o, the component of a spherical tensor of rank 2. Coherences ± 4 instead, grow slowly with r 
and dominate only within restricted intervals. 

At a preparation time near 60/xs the ratio between coherences 2 and changes drastically and their amplitudes 
become comparable to the rising 4-quantum coherence. Notice that this is precisely the preparation time at which 
C2(t2,£), the signal associated with the second QI, attains its maximum, showing that these states involve correlation 
between more-than-two spins. This characteristic indicates that these are multi- spin- correlated states. Notice that 
also near 110/is all the coherences attain comparable amplitudes, and this preparation time coincides with the third 
QI, -F 3 . 

We observed that the relative amplitudes of the various MQc components depend on which window along the 
acqusition time is used for calculating the MQc spectrum. This feature is shown in FigJTJb) where the time interval 
(along the acquisition period) at which the maximum value of each coherence component occurs is plotted as a 
function of the preparation time. The salient feature is that different coherences have a very different behaviour. 
The maximum contribution to coherences and ±2 on the X-basis come from the window near t m ax = 30 /is, except 
within the narrow preparation intervals which correspond to the second and third quasiinvariants. On the contrary, 
the maximum amplitude of coherences ±4 (on the X-basis) falls on the line t max = r, in other words, it behaves like 
an echo. This distinct behaviour agrees with a view that multi-spin correlations arise at longer times than two-spin 
correlations, and also that the echo like signals prepared with r > 70/is correspond to states of multi-spin nature. 

IV. DISCUSSION AND CONCLUSIONS 

In this work we present an experimental survey on the occurrence of several QI created using the JB sequence, in 
different dipolar-coupled spin systems. To account for the dependence of the JB signals on the preparation time, we 
propose a generalization of the method presented in refs.pjl. l37j. We demonstrate that multiple QI can be prepared 
in the dipole interacting proton spin system on LC molecules. The well-known dipolar order state, which is a two- 
spin object, comes at short preparation times, while the other QI which emerge sequentially at longer preparation 
times involve multiple-spin interactions. We propose a criterion to experimentally isolate each QI and oberve that 
the first three states in 5CB have evidently distinct spin-lattice relaxation times which implies that they are true 
quasiinvariants . 

Though at present the tensor structure of these quasi-invariants is still unknown, the multiple quantum coherence 
experiments demonstrate that the new states have multi-spin character, which implies that their tensor structure 
must involve products of individual angular momentum operators of many spins. Also, we verified that the observed 
echo-like behavior of the dipolar signal in LC for long preparation times is also a consequence of the occurrenece of 
many QI. 

In this way we demonstrate that the initial Zeeman order (of thermal equilibrium) can in fact be transferred to at 
least three dipolar constants of motion in LC (within the scanned short and intermediate r range), while only two 
constants of motion were observed in gypsum when B \\ [010] and only one at an orientation that makes the protons 
at the water molecules equivalent. The fact that compounds having very different numbers of interacting spins (11 to 
25 in LC and infinite in gypsum) admit the preparation of more than one dipolar QI, indicates that their occurrence 
is not exclusively determined by the size of the dipolar network. On the contrary, both the occurrence of multiple QI 
and the many-spin character they show, do depend strongly on the topology of the spin network, since it determines 
the quantum dynamics of the spin system. This assertion also agrees with the fact that the quality of our three QI 
model description varies from one LC to another. 

It is worth to note that the feasibility of transferring the Zeeman order to multi-spin quasi-equilibrium states 
is not an exclusive feature of spin clusters as the proton system in LC's as illustrated in this work. In fact, the 
"interpair" state observed in hydrated crystals admits multiple quantum coherences up to the fourth order on the 
X— basis, implying four-spin correlations 15]. Then the question arises about which are the necessary characteristics 
the dipolar network must have in order to admit the preparation of more-than-one dipolar QI. A coarse criterion 
could be the feasibility of truncating the weaker term of the dipolar Hamiltonian with respect to the stronger term, 
which generally implies the possibility of clasifying the dipolar couplings according to their intensity into "strong" 
and "weak" . In practice, this characteristic is also associated to the occurrence of a doublet in the NMR spectrum. 

According to our results, in LC's two-spin correlations dominate the coherent dynamics during the early timescale 
of the preparation period, however, there are narrow time windows where these correlations vanish and higher or- 
der correlations can efficiently give place to multi-spin QE states. Finally, for longer r values only the multi-spin 
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correlations are responsible for the echo-like behavior of the NMR signal. 
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FIG. 1: Stackplots of the JB signals as a function of the preparation time r in (a) Adamantane, (b) Gypsum and (c) MBBA 
liquid crystal. 
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FIG. 2: Acquisition time tmax at which the signals attain their maximum amplitude as a function of the preparation time r 
in the nematic phase of 5CB (top), 6CB (middle) and 8CB (bottom). Solid circles represent the experimental data. Open 
triangles and open circles correspond to t ma x(r) simulated from Eq.© with two and three quasiinvariants respectively. QIs 
selected as indicated in Table 1. 
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FIG. 3: Acquisition time t ma x at which the signals attain their maximum amplitude as a function of the preparation time r in 
a gypsum single crystal (top) and in nematic MBBA (bottom). Solid circles represent the experimental data. Open triangles 
and open circles correspond to tmaxir) simulated from Eq.Q with two and three quasiinvariants respectively. QIs selected 
as indicated in Table 1. Solid circles on the gypsum graph correspond an orientation B || [010], while open squares to an 
orientation that makes the protons at the water molecules equivalent. 
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FIG. 4: Signals Ti (dotted), Ti (solid) and Tz (dashed) for 5CB at 302K. 




FIG. 5: Detail of the experimental t max (r) in nematic MBBA (solid circles) and its description with Eq.Q considering two QI 
(open triangles) and three QI (open circles). 
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FIG. 6: Spin-lattice relaxation curves of three QI in 5CB at 302K prepared with different preparation times: r = 30/is, circles; 
r = 70/is, open squares and r = 134^is, full squares. 




FIG. 7: a) Maximum amplitude attained by each coherence, b) observation-time window at which the maximum amplitudes 
occur. 



